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I. INTRODUCTION 

A fascinating idea proposed by Wheeler in the early years of Quantum Gravity, is that, at the 
Planck scale, geometry may be bumpy due to quantum fluctuations. This is the quantum foam [T]. 
While intuitively natural, this idea is very complicated to put into action. In the present paper, we 
will use the framework of Quantum Graphity [SHE] to construct a simple model of quantum foam. 

A key feature of a quantum foam is its non-local nature. While non-locality is undesirable in 
quantum field theory, the situation in quantum gravity is open. It is often said that the only way 
to cure the divergences appearing perturbatively in quantizations of gravity without introducing 
new physics (i.e., string theory or super-symmetric extensions of gravity), is to introduce some kind 
of non-locality in the action which smears out Green functions evaluated on one point only. Until 
now, ghosts in the theory have blocked research in this direction (some progress has been achieved 
recently in [7]). For the purposes of the present work, it is important to note that there are two 
possible types of non-locality which contribute in different ways. One, violation of microlocality, 
disappears when the cut-off is taken to zero, while the other, violation of macrolocality, or disordered 
locality, does not [8J. Violations of macrolocality amount to the presence of what a relativist 
would call a wormhole [2], a path through spacetime disallowed in a Lorentzian topology. General 
relativity allows for such paths and, in principle, they should be taken into account in a full quantum 
theory of gravity. In principle, in order to have traversable wormholes, the common positive-energy 
conditions and some other conditions on the throats have to be satisfied. 

On the other hand, in graph-based quantum gravity states, such as in Loop Quantum Gravity 
[9], Causets |10j or Quantum Graphity [llj . spacetimes which are not macrolocal are very natural, 
and violation of macrolocality appears in the form of non-local links. A first study of the physics 
of these non-local links was carried out in [HI I12j. 

We propose, in the present paper, to use the framework of Quantum Graphity to provide a 
concrete implementation of Wheeler's quantum foam, based on the assumption that the non-local 
link can be used to cross from one end to the other one. 

Quantum Graphity models [3 , 4 are spin system toy models for emergent geometry and gravity. 
They are based on quantum, dynamical graphs whose adjacency is dynamical: their edges can be 
on (connected), off (disconnected), or in a superposition of on and off. We can interpret the graph 
as pregeometry (the connectivity of the graph tells us who is neighbouring whom). A particular 
graphity model is given by such graph states evolving under a local Ising-type Hamiltonian. The 
graphity model of [I], for example, is a toy model for interacting matter and geometry, a Bose- 
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Hubbard model where the interactions are quantum variables. 

In [5] , we solved the model of [I] in the limit of no backreaction of the matter on the lattice, and 
for states with certain symmetries that are natural for our problem, which we called rotationally 
invariant graphs. In this case, the problem reduces to an one-dimensional Hubbard model on a 
lattice with variable vertex degree and multiple edges between the same two vertices. The proba- 
bility density for the matter obeys a (discrete) differential equation closed in the classical regime. 
This is a wave equation in which the vertex degree is related to the local speed of propagation of 
probability. This allows an interpretation of the probability density of particles similar to what is 
usually considered in analogue gravity systems: matter inside this analogue system sees a curved 
spacetime. 

We will extend these results we obtained in order to describe a quantum foam: instead of a 
classical background state (a single graph), we will consider a state that is a superposition of many 
graphs. This amounts to a quantum foam with a superposition of Planck scale sized non-local 
links. In our setting, the intrinsic discreteness of the graph sets the minimum scale. Assuming 
foliability of the graph, we can define a metric distance as in [5]. We can then study the effect of 
the quantumness of the graph on the dispersion relations. 

Quantum Graphity models are lattice models in which the lattice becomes a quantum object. 
As in any lattice model, the continuum limit is obtained as in any other lattice theory, but consider 
it together for all the states on the graph. 

It is natural to construct graph states in which the largest contribution comes from the graph 
with the Lorentz invariant dispersion relations. The states with non-local links violate macrolocality 
and give corrections to Lorentz invariance. We will construct states with a distribution of non- 
local links which is suppressed by their combinatorial length. These states resemble coherent 
states as considered in Loop Quantum Gravity. In principle, they could be obtained as correction 
to the ground state due to a non-zero temperature bath in Quantum Graphity. The distribution 
depends on their density. We will then calculate the effect on the Lorentz dispersion relations in 
the continuum limit. The result is, as expected, a non-local differential equation for the evolution 
of the particle probability density. 

It is reasonable to expect that a non-local link will violate local Lorentz invariance. A particle 
can hop through the link and behave like a superluminal particle. As we will see, the presence of all 
these shortcuts has an effect on the local speed of propagation of probability density. Also, we will 
find that the probability density acquires a mass which depends on the density of non-local links. 
The overall dispersion relation is thus Lorentz invariant and with a square-positive mass. However, 
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this depends on the distribution and thus we will study two particular cases. Using the framework 
of Quantum Graphity and the techniques developed in [5], we will calculate the emergent mass 
and the constants appearing in the effective equation. 

This paper is organized as follows. In section II, we summarize the Quantum Graphity frame- 
work and the results of [5]. In section III, we show the effect of a superposition of graphs on the 
differential equation governing the time-evolution of the probability density. In section IV, we 
introduce our choice of the quantum state of the graph. In section V, we analyze two particular 
non-local link distributions and their effect on the dispersion relations. Conclusions follow. 



In the following we review the model, as defined and first studied in [4] and the effective geometry 
encoded in the graph, as obtained in [S]. 



In this section we will introduce briefly the model. For more a more detailed introduction we 
refer to the previous papers [H [5] . 

We associate a Hilbert space !Hj to the degrees of freedom on the nodes a graph, with i labelling 
the nodes. These degrees of freedom represent matter on the graph and thus can be, in principle, 
generalized to other fields. We choose !Hj to be the Hilbert space of a harmonic oscillator. We 
denote its creation and destruction operators by b\ and hi respectively, satisfying the usual bosonic 
commutators. Our N v physical systems then are N v bosonic modes and the total Hilbert space of 
such modes is given by 



II. THE MODEL 



A. Bose-Hubbard model on a dynamical lattice 



tosons — 



(1) 



i=l 



If the harmonic oscillators are not interacting, the total Hamiltonian is trivial: 




(2) 



8=1 



i=l 



The Hamiltonian reads as 




(3) 
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where h e is a Hermitian operator on Hi Hj representing the interaction between the system i 
and the system j. 

We introduce a primitive notion of geometry through the adjacent matrix A, the N v x N v 
symmetric matrix defined as 



1 if i and j are adjacent 

(4) 

otherwise. 



A defines a graph on iV„ nodes, with an edge between nodes i and j for every 1 entry in the matrix. 
The total Hilbert space for the graph edges is then 

JV„(JV„-l)/2 

^graph = (^) ^W-ej (5) 
e=l 

with JC e = Span{\0), |1)} a qubit representing on/ off links. Therefore, the total Hilbert space of 
the model is 

•W- — ^-bosons ® ^graphi (6) 

and a basis state in H has the form 

= \-q,(bosons)j ^ ^(graph)^ ^ 

= \m, ...,n Nv ) \ei, ...,e Nv ( Nv -i) ). (8) 

2 

The first factor tells us how many bosons there are at every site i, while the second factor tells us 
which pairs interact. 

We note that it is the dynamics of the particles described by 

^ho P = -E hop ^2Ai j (a\a i + /i.e.), 

i<j 

that gives to the degree of freedom \e) the meaning of geometry and h.c. denotes the hermitean 
conjugate. 

The hopping amplitude is given by t, and therefore all the bosons have the same speed. Note 
that, for a larger Hilbert space on the links, we can have different speeds for the bosons. 

As mentioned above, the long-term ambition of these models is to find a quantum Hamiltonian 
that is a spin system analogue of gravity. In this spirit, matter- geometry interaction is desirable as 
it is a central feature of general relativity. The above dynamics can be considered as a very simple 
first step in that direction. 



FIG. 1: A planar graph which is rotational invariant. 



In the present work, we study the model for a particular class of graphs that have been conjec- 
tured to be analogues of trapped surfaces. We are interested in the approximation k <C t, which 
can be seen as the equivalent of ignoring the backreaction of the matter on the geometry. As in 
|5], we will consider an Hamiltonian of the form 

H = H v + iThop- (9) 

In this case, the total number of particles on the graph is a conserved charge. H v and Hn n k s are 
constants on fixed graphs with fixed number of particles. The Hamiltonian is the ordinary Bose- 
Hubbard model on a fixed graph, but that graph can be unusual, with sites of varying connectivity 
and with more than one edge connecting two sites. Our aim will then be to study the non-local 
and quantum corrections to the effective geometry which can be encoded in the graph, as shown 
in [5]. Even on a fixed lattice, the Hubbard model is difficult to analyze, with few results in higher 
dimensions. It would seem that our problem, propagation on a lattice with connectivity which 
varies from site to site is also very difficult. Fortunately, it turns out that for our purposes it is 
sufficient to restrict attention to lattices with certain symmetries and then to restrict to an effective 
1+1 dimensional model. 



B. Rotationally invariant graphs and the encoded geometry 

Let us present next our definition of rotationally invariant graphs, which allows to reduce the 
problem to a 1-dimensional Bose-Hubbard model in the single particle sector. 

A graph G is called N -rotationally invariant if there exists an embedding of G to the plane that 
is invariant by rotations of an angle 2tt/N. In principle, the edges of the graph, once embedded, 
can be overlapping. The main property of the rotationally invariant graphs is that groups of sub- 
graphs can be labelled by an integer number i. These graphs can be very far from triangulations, 
as the rotationally invariant graphs in Fig. [T] and [2] show. 

These graphs can be labelled by a set of two integer coordinates, (n, 6 n ), where n labels a set of 
nodes, while 9 n is a coordinate internal to the subgraph. For convenience we will drop, since now 
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FIG. 2: A non-planar graph which is rotational invariant. 

on, the sub-index n in the 9 coordinates. 

We can make use of the coordinates (re, 9) in order to write the Hamiltonian defined by a 
rotationally invariant graph as 

N-l 

H rot = A nn ,b\ b n i e + h.c. 

0=0 n,n' 
N-l N-l 

0=0 <p=l n,n' 



(10) 



where b\ e {b n ,e) is the creation (annihilation) operator at the vertex (n,9), A nn i is the adjacency 
matrix of the graph and is the adjacency matrix of two angular sectors at an angular distance 
ip in units of 2n/N. 

Let us introduce the rotation operator M defined by 



Mb nfi = b r) e+iM 
Mbl e = b\ e+1 M. 



(11) 



The effect of the operator M is particularly easy to understand in the single particle case: 

M\n, 9) = Mbi fi \0) = bl fi+1 M\0) = \n,9 + l), (12) 

where we have assumed that the vacuum is invariant under a rotation M|0) = |0). Note that M is 
unitary and its application N times gives the identity, M N = 1. This implies that its eigenvalues 
are integer multiples of 2ir/M. 

Another interesting property of M is that commutes both with the rotationally invariant Hamil- 
tonians and with the number operator N p , 



[H rot , M] = [N p , M] = [H Iot , N p }=0. 



(13) 



Therefore H rot , N p , and M form a complete set of commuting observables and the Hamiltonian 
is diagonal in blocks of constant M and N p . In this sector of the Hamiltonian, we can reduce the 



s 



Hamiltonian to: 

L-l 

Ho = Y^ fn,n+l (|n)(n + 1| + \n + l)(n|) + ^^ n |n)(n| , (14) 

n=0 n 

with / nj „,+i depending on the degree of the graph and n being the label of the shells we are reducing 
with the rotational symmetry and L the total size of the one-dimensional lattice. 



C. Restriction of the time-dependent Schrodinger equation to the set of classical states 

Since we want to study the dynamics of a single particle on a fixed graph, it is only necessary to 
consider the single particle sector. The one dimensional Bose-Hubbard model for a single particle 



reads as in (14), where f n ,n+i are the tunneling coefficients between sites n and n + 1, p n is the 
chemical potential at the site n, and M is the size of the lattice. 

In this setup, let us introduce the convex set of classical states -Mc> parameterized as 

L-l 

p(t)=p(*(t)) =Y^^n\n)(n\, (15) 

n=0 

where ty n is the probability of finding the particle at the site n. The states in Aic are classical 
because the uncertainty in the position is classical, that is, they represent a particle with an 
unknown but well-defined position. 

Since our particle is under the effect of a noisy environment, its density matrix is going to be 
constantly dephased by the interaction between the particle and its reservoir. For a more detailed 
discussion about this procedure and the connection with the physics of decoherence, we refer to 
[5]. The dephased state in the position eigenbasis that best approximates p(t + At) can be easily 
determined by computing the double commutator of the previous equation, which was shown to 
lead to a closed equation in [5]. It obeys the evolution 

y£f *„(*) = £-l,n (*n+i(i) + *n-l(*) " 2¥ fl (t)) 
+ Cf„ 2 +l,n " /n-l,n) (*n+l(*) " *«(*)) • 

This equation becomes a wave equation in the continuum, 

d^(x,t)-d x (c 2 (x)d x v(x,t)) =o, (16) 

where 



h 2 



c(x) V 2p(x)El p £ hopV /27^ 



(17) 
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and ty(x,t) and f{x) are the continuous limit functions of ^> n {t) and f n ,n-\ respectively. Eqn 



( 16 ) is the equation of motion for a scalar field with a space-dependent refraction index. As it 
is well known, this equation in higher dimension is connected with the Gordon metric. In fact, 
to the refraction index it is possible to encode a space-time geometry with spatial curvature and 
no extrinsic curvature, i.e. a preferred direction of time. The time direction is the same of the 
quantum mechanical underlying model. This equation is the starting point for what we will do in 
the following. However, let us first recall how the continuum limit is performed. 

D. Dispersion relation and continuum limit 

Let us consider in more detail the translationally invariant case in which / n _i n = / and fj, n = fj, 



for all n. In this case, the continuous wave equation (16) becomes 

d^(x,t)-c 2 d 2 x O>(x,t) = 0, (18) 

where c is the speed of propagation. 

Let us recall how this limit was performed in [5]. Let us first introduce a discrete Fourier 
transform in the spatial coordinate and a continuous Fourier transform in the temporal coordinate, 
given by 

1 L_1 

¥„(*) = —^ k (t)e- 1 ^, (19) 

and *&k(t) = Ae luJkt + Be~ luJkt . After a straightforward calculation, we find that the relation 
between u) k and k is given by 



V L 

Now we can rescale uj k -4 0Jk/L (or equivalently c) and find that 



u k c = y/2 i/l- cos \^k). (20) 



u k c = LV2J 1 - cos H^-k ) . (21) 



and, therefore, 



lim £jk(L) ~ 2-7T-. 

L— >oo C 

That is, only modes that are slow with respect to the time scale set by c see the continuum. Note 
that by rescaling the speed of propagation c, the continuum limit can be obtained by a double 
scaling limit, -Ehop ~~ > Ehop/L and L — > oo for lattice size L. In this limit, the probability density 
has a Lorentz invariant dispersion relation. 
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III. A NON-LOCAL STATE DISTRIBUTION 



In this section, we show the effect of having a quantum superposition of graph in (24) on the 



equation (16). 



A. The effect of a quantum superposition of graphs 

In order to do the explicit calculation, we will modify the Bose-Hubbard interaction. Let us 
consider a one-dimensional Bose-Hubbard of the form, 

H = ^ j A i ^ 1 {a\a i + h.c.) (22) 

i 

and then consider its generalization, from Aij = <5j,i_i + to A{j = Nij, with JVy = Qjhj and 

bij,b\j the ladder operators on the Hilbert space of the link ij. is then the number operator on 
the Hilbert space of the graph, as usually considered in Quantum Graphity. This allows, instead 
of using fixed classical graphs, fixed quantum graphs, where the state \ipgraph) is superposition of 
different graphs. The full quantum hamiltonian for the system is, as usual, on an Hilbert space of 
the form 

H total) = Span{\il) graph ) ® \ip bosons )}. 

Using this, we want now to repeat the same calculation we performed in the previous paper, 
i.e. compute: 

dfMt) = ~i Tr{[H, [H,p(t)]]N' z }, (23) 

with ip n = (N^), N' n number operator on the bosons defined on the node n, and p the density 
matrix on the total system. Let us assume that the graph is not dynamical. We will also to use 
the Born approximation, that is, 

P(t)^p g ®p b (t), (24) 

with p g the density matrix of the graph and pb(t) the density matrix of the bosons. This ap- 
proximation allows us to consider a particle disentangled enough from the graph to be "followed" 



using the equation (16). It is also a physical requirement, which accounts for the existence of the 
particle on its own. In general, we expect that at long times the full hamiltonian thermalizes to a 
specific graph, depending on the parameter of the Hamiltonian which defines the metastable state. 
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Later on, we will rescale the coupling constant of the hopping Hamiltonian in order to obtain the 
continuum limit. Thus, one could think that this rescaling affects the state of the graph at infinity. 
However, the hopping of the bosons allows the graph to thermalize, as it has been shown in [3]. 
Rescaling this constant, just changes the time it takes for the system to thermalize, but not the 
asymptotic state of the graph. As a matter of fact, we do not know yet a Hamiltonian which gives a 
specific graph state asymptotically. However, the results of [13] in two dimensions and those of [3J, 
support the conjecture that, in general, such a Hamiltonian exists. For the time being, it is fair to 
say that the ground state of Quantum Graphity coupled to a thermal bath are rotational invariant 
graphs |6J. Thus, these graphs can at least be generated by a known effectively 2d-dimensional 
model. 

Based on these considerations, we conjecture the following graph state, \ip gra ph) = H>d) + IVw) 
with (tpci\4'ni) = 0. \ipni) is a correction to the classical graph state \ip c i) considered in [5] that we 
will discuss (and construct) in the next section. For the time being, let us consider the effect of 



this correction on eqn. (16). We have p g = \ip gra ph) (ip graph] ■ Thus: 



Pg = \lpcl}{lpcl\ + \tpnl){tpnl\ + (\lpnl) {4>cl\ + \lpcl) (25) 

Let us now evaluate these traces. A straightforward calculation shows that, 
E 2 

- -Srdfy n = Tr {{H 2 p + pH 2 - 2HpH)N z } 
n 

= 2 ^2 [Tr{AijAm n pg}Tr{alajdh^PbNz} 

ij,mn 

- Tr{A ij p g A mn }Tr{a} i a j p b a\ n a n N z }]. (26) 

We now substitute the equation for p g , and obtain: 

E 2 

- J ^ L dU n = Ki, n (t) + C n (t), 
n 

with Aip n (t) is the discrete second derivative and C n (t) is: 

C n (t) = 2Y, ijtmn [PijmnTr{a\aja\na n p h N z } 

-QijQmnTr{a\ajp b alna n N z }] , (27) 

with: 

Qij = (lpnl\Aj\tpnl), 

where we used the orthogonality condition (ip n i\ip c i) = 0. 

Our task now is to evaluate these two quantities on different classes of interesting states. 
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FIG. 3: The intuitive picture of non-local links inserted in the graph. 

IV. THE CHOICE OF THE QUANTUM STATE FOR THE GRAPH 

Let us now introduce the states on which we will evaluate the quantities defined in the previous 
section, Pij mn and Qij. Motivated by the fact that we can reduce using the translationally sym- 
metric graphs to one line, we will restrict our attention to a one-dimensional lattices. These states 
will resemble coherent states as considered in Loop Quantum Gravity. In principle, they could 
be obtained as correction to the ground state due to a non-zero temperature bath in Quantum 
Graphity. 

Let us consider first a metric on the classical graph, with d(i,j) the distance between the nodes 
of the classical graph \ip c i), with all the ordinary properties of distances. On a one- dimensional line 
this distance could be, for instance, given by \i — j\. Let us then construct states with non-local 
links on top. We want to penalize states with too long non-local links. We then introduce a factor 
p(i,j), which depends on a distance d(i,j) evaluated on the base graph, assuming that d(i,j) > 0, 
and a parameter I describing how non-local the links are w.r.t. the length of the graph. Then we 
define the operator: 

f l = Y,p(iJ)4p (28) 

i<j 

with 

£>(i 5 j) 2 = l, (29) 

i<j 

which ensures that p(i,j) 2 can be interpreted as a classical probability distribution. When applied 
to \ipd) this operator generates a superposition of all the possible non-local links which can be 
created on \vp c i), with a factor that with the distance of the links, 

IV4> = %d>, (30) 
and we can imagine to apply this operator several times to create more non-local links, 

\^) = %m. (3i) 
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The meaning to give to / is thus that of a cut-off in the length of these non-local links. Note that 
we can bias the number of links on which we want to peak the quantum non-local state the same 
way, 



v K s 



7? = E \ f i = eKTl ~ L ( 32 ) 

8=1 S - 

We see then that we can write the quantum state for the graph in the convenient form: 

Wnl)=[l + e Kfl ]Wcl). (33) 

This state depends explicitly on two parameters, I and K, and on the classical graph together with 
its distance. On this state we now want to evaluate: 

Pijmn = (lpnl\AjA mn \i) nl ) = (^ d \7f ^ A^A m ^jf \lp d ) , (34) 

and 

Qij = ^nl\Aij\^ n l) = Wdfi^AifiFWcl)- (35) 

Let us then consider first the average. We note that, since Aij acts like a projector, and states 
with different powers of the 7} operators are orthogonal, we can write: 

Q*i = E -^WcllTtAjTf \^ d ). (36) 

s=l 

To clarify the idea, let us consider the case in which we add just a link. In this case, the 
state is the sum over all possible links which can be created, with a factor p 2 (i,j).This link can 
be created in one way only, and so Aij projects on the only state which can be non-zero. A very 
straightforward calculation shows that 

{4> cl \T? 1 A ij T l 1 W cl ) = 2p 2 (i,j). (37) 

For the higher order term, we instead have: 

s-l 

^ cl \T} s A ij T l s \^ d )=p\i,j) E UAii,3i)- (38) 

h,ji,— ,is-i,js-i 1=1 



It is easy to see that 



E l[p 2 (ii,ji)*2 s s(lL) s , (39) 
h,ji,--- ,i s -i,js-i 1=1 
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due to the fact that the integration is over the line, while the distribution has an extension of circa 
I combinatorial points. The factor 2 s comes from the fact that there are 2 points we are summing 
over and the s factor from the s sums appearing in Tf. Thus, we can write: 

(falTl'AijTflfa) = c s p 2 (i,j) 2 s s (L (40) 



In principle, given a distribution, we can calculate this factor from eq. (39). We will calculate 



these factors later for two particular distributions. Plugging eqn. (40) into Qy, we obtain 



00 jy;2s 

Q*i = E^T 25 s ( L l) s - 1 p 2 (hj)cs = p 2 (i,j)R(K,l L), (41) 



with: 



00 f(2s 

R(K,l L) = Y,^2~ c s 2 s s (I L) s ~\ (42) 

and, therefore, 

QijQmn = p 2 (i,j)p 2 (m,n)R(K,l L) 2 . (43) 
We can, in fact, do an analogous calculation for Pij mn and find that: 

Pijmn = p 2 {i,j)p 2 (m,n)L( y K,l L), (44) 

with: 

L(K, IL) = J2 -|2-C(Z L Y~ 2 2 s s. (45) 

8=1 S - 

Going back to the original problem, we find that the correction to the discrete Lorentz equation 

is: 

Cz = 2Y, ij>mn p 2 (i,j)p 2 (m,n)[L(K,l L)Tr{a\afcLhnP b N z } 

-R(K, I L) 2 Tr{4 aj p b ala n N z }] . (46) 

If we define: 

S(K, I L) = (I L) R(K, I L) = (I L) 2 L(K, I L), (47) 
with S(K, I L) = YZLx Cs L Y 2" s > then we obtain: 
C z {t) = 2 Y, p 2 {hj)p 2 (m,n)S{K^L)[Tr{a\ aj ala n p b ^ (48) 

ij,mn 
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FIG. 4: A plot of the function S(x) which appears in eq. (48) 



We see that the function S(K, I L) depends, as a matter of fact, on £ = K\Jl L, S(K, I L) = S(£) = 
Yl^Li °s [§r] 2 2 s s. A plot of the function S(K,IL) can be found in Fig. |4j The traces can be 
evaluated, as done in ([5]), and the result is: 



C z (t) = ^p 4 (^, L)ty z - S(K,l L)ipj). 



(49) 



Some comments are now in order. First of all, note that the equation has the shape of a second 
derivative. To understand this, we can look at a term of the form ^|fc|>2 J(k)(i>z ~ i^z+k)- This 
term can be written as: 

J(k) ■■■ = - J(2) (V>,+2 - 2V* + i> g - 2 ) - J(3) (^+3 - 2i/} g + ^_ 3 ) - J(4) • • • . 



This is a sum of discrete second derivatives with a non-local mass, 

fc-i 



J (k) (tpz+k - tyz + ipz-k) = -J(k) {i) z+ k~i + 4>z-k+i) -J(k)^2 (tpz+i+i ~ 2ipz+i + 4>z+i-l) , 

i=2 

so we expect, in the end, to obtain a mass term out of this equation and, when we will have 
rearranged all the terms, we will. 

Note that, for the case c s = 1, = 1 for £ = 0.903, and so K = ^S. We then see that K 2 
plays the role of the density of non-local links per units oil L. 

To end this section, we have to calculate the norm of this state. This can be written as: 



\^nl\i>nl)\ = ^l + (i>d\e KT leKT l ^ d )2 + 2U{{i , d \ e KT l \ i , cl)h 



(50) 



which reads, 



\^Mnl)\ = V 1 + {&d\e KT > e^ d )2 _ i) } 



(51) 
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and substituting for the 7 operators, we finally find 



^ 1 + ( E pa E ^(^-) 2 ) = y 1 + ( E (^2 2Ss ) = + 



(52) 

We can thus normalize the graph state by dividing by a factor of N. 

V. THE MODIFIED DISPERSION RELATION DUE TO DISORDERED LOCALITY 

The general case. We will now discuss the continuum limit. As we have seen, the continuum 
limit is obtained by rescaling Eh op — > Eh op /L and then sending L — > oo. Please note that Eh op 
appears whenever we hop with a particle, so in these calculations it appears everywhere but in 
the df term. In order to perform the continuum limit, first we have to make sense of the quantity 
(ip z — S(K,l L)ipj) at least for the flat case, which we know correspond to Lorentz from [5]. We 
can add and subtract, 

(S(K, I L) - 1)^ + (S(K, I L)^ z - S{K, I L)ijjj) = 
= (S(K, I L) - 1)^ + S(K, I L)(ip z - + + </>,■). 

In the continuum limit this becomes (S(K,l L) — l)if)(z,t) + S(K,l L) J* ; d%rf}(£, t)d£, and thus 
C z (t) reads: 

f r 4/ sr (S(K,l L)-l)S(K,l L) N S 2 (K,l L) f* 
C z (t) = J^dx p\z, x ){±-^ 1 ^ V + ^ ' ' J d^(tm\, (53) 



which is: 



4 , ,(S(K,l L ) - l)S(K,l L) S 2 {K,l L) f 4/ 
WW + ~W 2 j L " ' 



C z (t) =4>(z,t) I dx p (z,x) -rj-jr^ 1 n ^ 2 J^p (z,x) J d£ dx. 

(54) 



This can be written as: 



C z (t) = xp{z, t)F(K, I L) + 0(K, I L) / p A (z, x) / d&(£, t) d£ dx, 

with F(K,l L) = J L dxp 4 (z, x) {S{K ' 1 L )~ljf {K ' 1 L) ,0{K, I L) = ■ Please note here that 

these steps have been performed naively, though we have an explicit dependence on L in S. It is 
important to point out that the only way to keep the function S(l L, K) finite is to rescale the 
quantity K 2 l ~ ^t~- To keep the discussion simple, let us discuss this point at the end of the 
section. L is the combinatorial length of the 1-d lattice we are considering, and over which ip(x,t) 
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is defined. Thus the equation of motion for the flat case is given, in the continuum, by: 

[d 2 - c 2 (l + S 2 (K,l L))d 2 z - F(K,l L)]rp(z,t) = 0(K,l L) J p 4 (z,x) J d&fat) d£ dx, 

which is an integro-differential equation for the field integrated over the line, which shows the 
strong non-local character of the equation. 

We note that there is a contribution to the speed of propagation of the signal, due to the fact 
that particle can hop on many more graphs than the single classical one. This factor contributes 
with a c 2 S 2 (K,l L) added to the effective speed c 2 . Let us stress that this contribution is merely 
due to the fact that there are many more graphs in the superposition, and not due to the fact that 
the particle can hop further: this is kept track of in the C z (t) term of the equation. Also, we see 
that F(l,K) becomes a mass, due to non-locality, while on the r.h.s. there a new term appears. 
We can further reduce the equation by evaluating the integrals. It is clear that in order to have a 
finite result, which is physically expected, we have to rescale at this point only I ~ l/L, keeping K 2 
independent from L. Said this, we see that the distribution itself, when is well chosen, becomes a 
5 function and therefore the models becomes local again. 

Let us now calculate the terms at the leading order in l/L, since that is what we are interested 
in. The discrete differential equation becomes: 

L 

[d f _ c 2(i + s 2 (K,l L))d 2 - F(K,l L)]^ z (t) = ~0(K,l L) ^ p A (z, x)tp z (t), (55) 

x=0 



where d 2 is the discrete spatial second derivative. Using now (19), we see that the dispersion 
relation for the field becomes: 



/ 2tt 

w fc c(l + S 2 (K, I L)) = V2J1 - cos (— k) + F(K, I L) + p A (k)0{K, I L) 

Please note that with this rescaling of K, we have that S(K, I L) can be expanded in even powers 
of l/L: 

Thus, we see that the superluminal effect, which is, the factor 1 + S 2 (K, I L), becomes one in the 
limit L — > 00; also, in the same limit, only the part quadratic in K survives. At this point the 
equation would become, in the continuum: 

[df - c 2 d 2 - F(K, I L)]i;{z, t) = -0{K, I L) j p 4 (z, x)t/j(z, t) dz (56) 
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with F(k,l L) = F(K,l L) + 0(K,l L) j L p 4 (z,x) dx. Note that, while F might seem to be 
dependent on the point z, being F dependent on z—x and integrated over x, it is indeed independent 
from it. In particular, if we define I L = £, in the limit L — > oo and with the rescaling of K and I, 
S(K, I L) — > 2K 2 £. We see now that the only way to obtain the continuum dispersion relation by 
rescaling c — > c/L, as done for the single-graph state, is to rescale also K, with K — > K/L. 

Just as an exercise, we can insert a trivial spatially-constant solution, which then becomes of 
the form d^{t) = R(K, I L)4>(t). where R(K, I L) = F(K, I L) + 2 0(K, I L) f L p 4 (z, x) dx. Note 
that this quantity is always positive, so constant solutions are stable. Let us try to find a generic 
solution, instead. Let's do it for the equation: 

[8 2 t - c 2 d 2 x + c 2 (#(x, t) = -c 2 J a(z, x)4{y, t) dy. (57) 

Since the equation is linear in the field ip, we can solve it by means of a Fourier transform. We then 
look at the dispersion relation for the function ip(x,t), with q and P generic functions. We can do 
it by Fourier transform. In this case, the integral on the right, being a convolution, becomes just 
the product of the Fourier transform of the single functions. Thus we have: 

-u? + k 2 c 2 + c 2 q = -c 2 a(k), 

and we have that: 

u = ±c\Jk 2 + q + a(k). 

Now, of course a(k) depends on the distribution of non-local links that we inserted in the 
wavefunction of the graph. 

Two specific distributions. Let us consider two specific cases: 

1 rr _(x~y) 2 

• pi{X — y) = 7T4 VI e 2(2 ; 

• P2{x - y) = V^~l e^^^r 1 . 

In these cases we find, using standard tables of Fourier transforms: 
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and thus, keeping track of all the factors, we obtain: 



and 



with Fi(K,l L) = a/2 2 (K,l L) and F 2 (K,l L) = 2 (K,l L), which can be calculated by 
evaluating f L pf(x — y)dx. We have that S 2 (K,l L) = AK 4 ^ 2 /L 4 and thus can be neglected with 
respect to 1. Also, since the c contribute with a factor of L 2 within the square root, also S 2 can 
be neglected, and it contributes only the mass term in the L. Now we note a nice property: both 
the two distributions go to for k — > oo, that is, at high energy the dispersion relations become 
Lorentz again. We see then that the total effect the one of having an effective scale-dependent 
mass, which runs from one mass to another one, in both cases: 



1 / fc 2;2 

m i( k ) = FtV S 2 {K,l L) + S{K,l L)e- — , (60) 



m 2 (k) = ^ L \js 2 {K,l L) + S(KJL)* w ^- ¥ . (61) 



The masses which are intertwined are given by 



mi(0) = ^S 2 {K,l L) + S(K,l L), mi(oo) = L) , (62) 



m 2 (0) = ±-J&{K,l L) + -S(K,l L), m 2 (oo) = S &> 1 L \ (63) 
£ L V vr I L 

This property, of intertwining two different masses between k = and A; = oo is shared by any 
function which is at least C 1 . It is remarkable, instead, that the mass at k = oo does not depend on 
the distribution we inserted at hand. In fact, any C r distribution will lead to a Fourier transformed 
distribution which goes to zero at k = oo as l/k r and thus tend to a finite value for the mass. 
Note that, if we send I — > 0, as required to have S finite, the dependence on the scale seems to 
disappear, leaving a Lorentz dispersion relation with a mass which depends on the function S. 
However, we have to remember that, in fact these Fourier distributions come from the discrete 
dispersion relation. There, the distributions depend on 2-irk/L. If we define kL = k, then we have 
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FIG. 5: The running of m i (k) for x = 0.1, I = 1 and L = 1. 




FIG. 6: The running of mj(fe) for a; = 0.1, I = 1 and L = oo. 



that the distributions cancel out the dependence on L, leaving exactly (62) and (63) but dependent 



on this new momentum k. Still, this mass depends on the distribution we have chosen through 
di{k = 0) and so it has a valuable effect. We plot the running of (I L)rrii{k) as a function of 
x = K 2 I L, for the case L = 1 in Fig. ^ and L = oo in Fig. 

We would like to point out that the appearence of a mass which is square-positive is rather 
surprising. The physical reason is that, before starting the calculation, we would have expected 
that the presence of these non-local links would have shown a superluminal effect due to the 
non-local links themselves. However, the effective speed of propagation is higher because of the 
superposition of the graphs and not the non-local links. Indeed, the non-local links contributed 
only in the mass, thus the term C z (t) additional to the differential equation we obtained. Besides, 
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this mass is square-positive, thus it is an effective mass and not a tachyonic one, which we would 
have expected from the presence of non-local links on physical grounds. The fact that it is square 
positive comes merely from the fact that the equation comes from a quantum mechanical average, 
and thus the terms appear squared. 



VI. CONCLUSIONS 



One of the most striking theoretical consequences of General Relativity is the existence of 
wormholes and black holes. While the second is currently investigated experimentally, less is 
known about the first. Here we discussed something which in principle is very similar, quantum 
states which violate macro-locality. Besides, it could be that the quantum state of the Universe 
is a superposition of spacetimes with non-local links. In the present paper we considered such a 
possibility in a toy model constructed using the framework Quantum Graphity. In order to do 
so, we had to extend the results of [5] to a case in which the quantum state of the background 
is a superposition of many graph states. The superposition of these graphs was chosen so that 
it is dominated by a graph on which, as we showed in earlier papers, the expectation values of 
number operators of the bosons hopping on it satisfy a closed equation for probability density in the 
classical regime, i.e. a wave equation. We extended the formula previously obtained and studied 
a particular case: graphs which violate micro- and macro- locality. As discussed, a violation of 
macro-locality can be interpreted, within the model, as the presence of spatial non-local links in 
the background spacetime. This is a concrete example of a quantum foam within the framework 
of Quantum Graphity (3j 0] . The graph state was chosen on the basis of what we know from low 
energy physics, which is that Lorentz invariance is satisfied up and above the Planck scale[14|. 
We also used a class of graphs introduced in [5], rotationally invariant graphs. By exploiting 
their symmetry, the problem can be reduced to a 1-dimensional one, i.e. Bose-Hubbard model on 
a line with specific couplings depending on the connectivity of the graph. We thus constructed 
the states that are corrections to the low-energy physics by assuming that the non-local links 
are suppressed by a length according to a certain distribution. The length is measured by a 
combinatorial distance based on the low energy graph and which defines the state. We studied for 
the cases d(x, y) = (x — y) 2 and d(x, y) = \x — y\. We found that, in the continuum limit, there is 
no superluminal effect on the low-energy physics, i.e. the speed of propagation is intact. However, 
there is an appearance of a mass dependence on the constants of the distribution and that can be 
calculated within the model. These masses are square-positive and thus do not violate the physics 
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of the restricted Lorentz group, i.e., are not tachyonic. A simple analysis showed that this mass 
runs with the energy scale and, in particular, runs to zero at high energy. It is interesting to ask 
whether a similar phenomenon happens for the other fields. This analysis suggests the possibility 
that a quantum foam could contribute to the mass of a quantum field. As suggested in [8] and |12J . 
the possibility of having non-local link states within Loop Quantum Gravity is very natural. Also, 
it has been suggested that these states could contribute to the dark energy puzzle. The results 
of the present paper suggests that, as in [12], the quantum foam contributes to the mass of fields 
hopping on such a superposition of spacetimes. We believe that such possibility needs to be further 
investigated. 
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